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ABSTRACT 



Or 



We explore potential impacts of nuclear burning on assisting an onset of the 
neutrino-driven explosions of core-collapse supernovae. By changing the neu- 
trino luminosity and its decay time to obtain parametric explosions in one-(lD) 
and two-dimensional (2D) models with or without a 13-isotope a network, we 
study how the inclusion of nuclear burning could affect the postbounce dynamic s 
for four progenitor models ; thre e for 15.0 s t ars of iLimongi fc Chieffil (120061 ). 
Woosley h Weaver! (119951). a nd IWooslev et all (120021 ) and one for an 11.2 M & 
star of IWooslev et al.l ( 120021 ). We find that the energy supply due to nuclear 
burning of infalling material behind the shock can energize the shock expansion 
especially for models that produce only marginal explosions in the absence of 
nuclear burning. These models enjoy the assistance from nuclear burning typ- 
ically in the following two ways, whether the shock front passes through the 
silicon-rich layer, or later it touches to the oxygen-rich layer. Depending on the 
neutrino luminosity and its decay time, the explosion energy increases up to a few 
times 10 50 erg for models with nuclear burning compared to the corresponding 
models without. The difference in the explosion energy is generally smaller in 
2D than in ID, because neutrino-driven convection and the SASI in 2D models 
enhance the neutrino heating efficiency, which makes the contribution of nuclear 
burning relatively smaller compared to ID models. We point out that these fea- 
tures are remarkable only for the Limongi-Chieffi progenitor both in ID and 2D, 
which possesses a massive oxygen layer with its inner-edge radius being smallest 
among the employed progenitors, so that the shock can touch to the rich fuel on 
a shorter timescale after bounce. Considering uncertainties in progenitor models, 
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our results indicate that nuclear burning should remain as one of the unignorable 
ingredients to foster the onset of neutrino-driven explosions. 

Subject headings: supernovae: general — neutrinos — hydrodynamics — nuclear 
reactions, nucleosynthesis, abundances 



1. Introduction 



Ever since the dawn of modern core-co llapse supernova (CCSN) theory, the neutrino- 
heating mechanism fjColgate fc White! Il966l ). in which a supernova shock is reviv ed by neu- 
trino energy deposition to trigger explosions ( 1Wilsonlll985t iBethe fc Wilsonlll985l ). has been 
the leading candidate for the explosion mechanism for more than four decades. However, 
the simplest, spherically-symmet ric (ID) form of thi s mechanism fails, except for super- 



AGB stars at the low-mass end flKitaura 



flRampp fc JankakoOO ; ILiebendorfer et al. 



et a 



2001 



2006 



, to explode canonical massive stars 



Thompson et all2003l ; ISumiyoshi et al.ll2005 ). 



Push e d by accumulating s upernova observations of the blast morphology (e.g. JWang fc Wheeler 



20081 ; iTanaka et al.l 12012k and references therein), a number of mult i- dimensional (multi- 



D) hydrodynamic simulations have been reported so far , which gives us a c onfidence that 
hydrodynamic motions associated with convection (e.g 



Herant et al. 1992; Burrows et al. 



19951 : IJanka fc Murlerill996l:lFrver et al 



Instability (SASI, e. g., 



2007 



Foglizzo et al. 



2009a 



Fernande 



Blondin et al 



2006 



2010 



2007. 



2012 



Hanke et al. 



2002 



2003 



Frverll2004j) and the Stand i ng-Accretion-Shock- 



IScheck et al 



Iwakami et al. 



20111 ; 



2004, 



2008 



Miiller et al. 



neutrino-driven explosion (see collective references in |Jankall2012l ; iKotake et al.ll2012l ). 



2006; 



2009 



Ohnishi et al. 



2006 



Fernandez fc Thompson 



201 2bl) can help the onset of the 



In fact, the neutrino-driven explosions have been obtained in the following first-principle 
two-(2D) and three-(3D) dimensional si mulations in which t he spectral neutrino transport is 
solved by various ap proximations (e.g., Kotake et al. ( 20121)). Among them are the work by 



the Garching group teuras et all I2006a1 b1: 



Marek fc Jankall2009l ; iMiiller et al.l 120 111 l2012al ) 



who included one of the best available neutrino transfer appro ximations by the ra y-by- 
ray variable Eddington factor method, by the Oak Ridge group ( iBruenn et al.ll2010l ) who 
included a ray-by-ray multi-group flux -limited diffusion transport with the best ava ilable 

who 



weak interactions, and by our group ( ISuwa et al.l l2010l 12011 



Takiwaki et al 



employed a ray-by-ray isotropic diffusion source approximation (ILiebendorfer et al. 



2012 



20091 ) 
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with a reduced set of weak interactional- 



This success, however, is accompanying further new questions. First of all, the explosion 
energies obtained in these 2D simulations are typically underpowered by one or two orders 
of magnitudes to explain the canonical supernova kinetic energy (~ 10 51 erg, see table 1 



in 



Kotakd (120111 ) for a sum mary). Moreover, the so ftest version of Lattimer & Swesty 



(LS) equation of state (EOS) (ILattimer fc Swestylll99ll ) with an in compressibility at nu clear 
density, K, of 180 MeV, have been often employed (see, however, iMiiller et al.l (j2012al )). A 
serious flaw is that the LS180 EOS cannot support a 2M^ cold neutron star which certainly 
exists in the universe ( Demorest et al. 2010 ; Steiner et al. 2010 )^1. 



What on earth is missing furthermore? 3D hydrodynamic s has been pointed ou t to 
foster the onset of neutrino-driven explosio ns compared to 2D (INordhaus et al.ll2010l ). al- 
though it is still under considerable debate ( IHanke et al.l 1201 ll ; iTakiwaki et al.l 120121 ). Very 
recently, ge neral relativity has b ee n repo rted to help the onset of multi-D neutrino-driven ex- 



plosion s by iMuller et a. 



.1 ( 2012al Jbl. I2OIII ) in 2D simulations with detailed neutrino transport 



and by iKuroda et al.l (120121 ) in 3D simulations but with approximate neu trino transport. 



Impac t s of nuclear EOSs in the multi-D con t ext h ave been investigated by iMarek fc Janka 
f|2009h : IMarek et al.l ffio09h and ISuwa et all (ESSE However, there may still remain fur- 
ther room to study more detailed nuclear physical impacts in these first principle multi- 
D simulations, s uch as the density dependence of symmet ry energy and the skewness of 



comp ressibility (jSteiner et al 



e.g. 



2010 



Sumiyoshi fc Ropkdl2008l ; lArcones et al. 



200 



Lattimer fc Liml 120121) and influence s of light nuclei 



Nakamura et a 



20091) and of inelastic 



neutrino-nucleus scattering (e.g., iHaxtonl Il988l ; lOhnishi et al.l 120071 ; lLanganke et al.l 120081 ) 
on enhancing the neutrino heating rates in the gain region. The neutrino-driven mechanism 
would be assisted b y other candidate mechanisms such as the ac oustic mechanism (e .g., 



Burrows et al 



20061) or the magnetohydrodynamic mechanism (e . g..lKotake et al.l (l2004alJbr); 



akiwaki et al. 



(l200il2009h:lrjjuTows et Jfco_07aj) 



(l201ll ): ITakiwaki fc Kotakd (1201 if ), see also 



Guilet et al. 



Kotake et al 



( |201ll ); IObergaulinger fc Janka 



J2006I ) for collective references 



therein). Other possibilities include QCD phase transition s in the core of the protoneu- 
tron star (e.g., iTakahara fc Satol Il988l; ISagert et al.l 120091 ) viscous heating by the mag- 



netorotational inst ability ( iThompson et al.l 120051 ) . or energy dissipation via Alfven waves 
JSuzuki et al.ll2008l ). 



1 See lSumivoshi fc Yamadal (|2012f ) for collective references about detailed neutrino transport schemes. 
2 The maximum mass for the LS180 EOS is about 1.8M (e.g.. lO'Connor fc Ottl |201ll ); lKiuchi fc Kotake 



(120081) ). 

3 see also ICouchl (|2012l ) in which the EOS impact has been studied in parametrized 2D models. 
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Joining in these efforts to look for some possible missing ingredients, we explore in this 
study potential impacts of nuclear b urning on fos t ering the onset of the neutrino-driven explo- 
sions. To the best of our knowledge, iJanka et all ( 120011 ) were the first to clearly point out that 
an additional energy released by nuclear burning of infalling material behind the shock could 
make a significant contribution to affect the explosion energy (see their Eq.(5)). The mass in 
the silicon (Si) layer, depe nding sensitively on the p r ogenitor masses and s tructures, is in the 
range of ~ 0.3 - O.6M JWooslev & Weaver! EjjSa IWooslev et all I2OO2I ; iLimonri & Chieffi 
20061 ). Since the release of nuclear energy per gram in Si burning is ~ 10 18 erg/g, the nuclear 
energy E m W , of a few 10 50 er g is expected to be deposited by the explosive nuclear burning. 



However, IJanka et al.l ( 1200 ll ) hypothesized that based on an order-of-magnitude estimation, 
the energy gain by the nuclear burning might not significantly affect the explosion dynamics 
for progenitors above 2OM , because the gravitational binding energy of the overlying, out- 
ward accelerated stellar layers outside the iron core, could be roughly the same order of E mc , 



possibly leading to the cancel lation of the two effects. The Oak Ridge group ( jBruenn et al. 



20061 ; iMezzacappa et al.ll2007l ) obtained explosions in their 2D radiation-hydrodynamic simu- 
lations for 11.2 M Q and 15.0 M Q stars, only when an alpha network calculation was included. 
They pointed out that the oxygen burning assists the (weak) shock to move farther out due 
to the additional pressure support in the vicinity of the shock. 

In the present work, we take the following strategy to investigate the possible influ- 
en ce of nuclear burn i ng on the neutrino- heat ing mechanism. Firstly we try, in the spirit 
of iBurrows fc Goshyl (119931 ) and iJankal (1200 ll ). to find a critical condition in ID, in which 
nuclear burning affects the criteria of explosion. For th e sake of our systemat ic survey, we 
employ a light-bulb scheme to trigger explosions (e.g., IJanka &: Miillerl Il996h . Previously 
the role of nuclear burning seems to be considered as negligible using a very limited set of 
progenitor models (for example, see discussions in lNordhaus et al.ll2010l ). but we will show 
that for a virgin progenitor model, nuclear burning can really push the weak shock farther 
out to help explosions. After we first give a detailed explanation in ID models how nuclear 
burning could affect the postbounce dynamics, we then move on to discuss our 2D models 
to study how nuclear burning interacts with multi-D hydrodynamics. 

This paper opens with the description of numerical setup including information about 
our hydrodynamic code with nuclear network and about initial models (Section 2). Results 
are given in section 3. We summarize our results and discuss their implications in Section 4. 
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2. Numerical Setup 
2.1. Hydrodynamics with Nuclear Network 



We solve the hydro dynamic equations corresponding to the conservation of mass, mo- 
mentum, and energy, 

P +pV-v = 0, (1) 



dt 



-Vp - pV$, 



de „ 
3i + V 



[(e + p)v] = -pv -V$ + p{H-C + Q), 



(2) 
(3) 



where p is the mass density, v the fluid velocity, p the pressure, $ the gravitational potential, 
and e the total energy density, respectively. The Lagrangian derivative is denoted by d/dt = 
d/dt + v ■ V. To treat Newtonian self-gravity, a monopole approximation is em ployed. The 
tabul ated realistic equation of state based on the rela tivistic mean field th eory (IShen et al. 
19981 ) is implemented according to the prescription in iKotake et al.l (120031 ). 



In this study, we employ the so-called light-bulb scheme (IJanka fc Mullerlll996l ). in which 
neutri n o hea ting and cooling is adjuste d parametrically to trigger explosions. Following 
Jankal (120011 ) and iNordhaus et al. feoioh . the neutrino heating (H) and cooling rates (C) 
are given by, 



H = IMA x 10 



20 



L„ 



C = 1.399 x 10 



x 



20 



10 52 ergs 1 

f \ — 2 



T 



100 km 
T \ ^ 



4MeV / 
{Y n + Y p ) e 



[erg/g s] 



[erg/g s] 



(4) 



(5) 



2MeV J 

where L Uc is the electron-neutrino luminosity that is assumed to be equal to the anti-electron 
neutrino luminosity (L Pc = L Uc ), T Vc is the electron neutrino temperature assumed to be kept 
constant as 4 MeV, r is the distance from the center, T is the local fluid temperature, Y n 
and Y p are the neutron and proton fract ions, and r v i s the electron neutrino optical depth 
that we estimate according to Eq. (7) in lHanke et al.l (120111 ). 



Note in this study th at neutrino luminosity is assumed to evolve exponentially with 
time (IKifonidis et al.l 120031 1 as 



L Uc = L Pc = L u0 exp(-t ph /t d ), 



(6) 
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where L u q denotes the initial luminosity, t p b is the time measured after core bounce td is 
the decay time, respectively. L u0 and td are treated as free parameters. Note that neutrino 
heating and cooling are switched on only after core bounce. 

We are keeping track of 13 species of a network (from 4 He to 56 Ni) by solving a separate 
advection equation for each species. Q in Eq. ([3]) denotes the net energy deposition rate by 
nuclear burning. The nuclear r eaction network is mainly based on t he REACLIB database 
(IRauscher fc Thielemannll2000l). E xperimentally determined masses (lAudi &: Wapstralll995l ) 
and reactions (lAngulo et al.l Il999l ) are adopted if available. It should be noted that our 
network does not include the photodissociation of iron elements because Shen EOS adopted 
in this study takes account of these endothermic effects. Note also that we solve the reaction 
network only for the grids where T < 5 x 10 9 K, assuming that above this temperature the 
local chemical composition is in nuclear statistical equilibrium (NSE). 

Only after core bounce, neutrino heating and cooling is switched on, according to the 



prescriptions (Equations (J4]) and (jSJ)) assu ming L 



and 7k = 7k. Before bounce, 



we employ the Y e prescription proposed by iLiebendorferl (120051 ). in which Y e is given simply 
as a function of den s ity, a n d after that, we r efrain from solving the c hange of Y e following 
Murphy fc Burrows! (120081 ) ; iNordhaus et al.l (120101 ) and lHanke et al.l (120111 ) (s ee, however. 



Ohnishi et al.l (120061 )). As for the hydro-solver, we employ the Z EUS-MP code ( Haves et al. 



20061 ) which has been modified for core-collapse simulations (e.g.. llwakami et al.ll2008l . 120091 ). 
The computational grid is comprised of 300 logarithmically spaced, radial zones to cover from 
the center up to 5000 km for our ID models with 128 polar (8) uniform mesh points for our 
fiducial 2D models. 

In order to induce non-spherical instability after the stall of the prompt bounce shock, 
we have added a radial velocity perturbation, 5v r (r, 9, 0), to the steady spherically symmetric 
flow according to the following equation, 

v r (r 1 6,<f ) ) = v r (r) + 6v r (e,<f>), (7) 

with 

5v r = 0.01 x rnum x v®(r, 9) , (8) 

where v®(r,6) is the unperturbed radial velocity and 5v r is the random multi-mode pertur- 
bation with a random number —1 < rnum < 1. 



2.2. Progenitor models 



In this study, we e mploy four progenito r mode ls; three for 1 5.0 stars of 

(l2006l . hereafter LC15). IWooslev fc Weaverl Jl995l . WW15), and lWooslev et al. 



Limongi Sz Chieffi 



, WHW15) 



(12002 
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Fig. 1. — Precollapse composition distributions for the 15_M© stars of iLimongi fc Chieffi 
fcflQd) (labeled by L C15, top left), IWooslev fc Weaver! f|l995h (WW 15, bottom lef t) , and 



Wooslev et all (120021 ) (WHW15, top right) and for the 11.2 M star of IWooslev et all (120021 ) 



(WHW11, bottom right). 
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and one for an 11.2 M & star of IWoosley et al.l (120021 WHW11). For all the models, Figured] 
shows the precollapse composition profiles near from the outer edge of the iron core to out- 
side. As we will explain in the next section, burning of the oxygen shell behind the (weakly 
propagating) shock plays an important role in assisting the shock expansion. Therefore, 
the earlier the oxygen layer touches the (stalling) shock after bounce, the better it could 
work. Having this in mind, let's see Figure 1 again. Among the three variants of the 15 M Q 
progenitors, the inner edge of the oxygen layer (seen as a sharp decline in solid red lines of 
Figure 1) is positioned much closer to the center for models LCI 5 (closest, top left panel) 
and WHW15 (next closest, top right panel) compared to model WW15 (bottom left panel). 
Table 1 shows a summary of the precollapse abundance distributions, in which each quantity 
from the left to right column corresponds to the different progenitor models, the progenitor 
mass, the mass of the iron core, the outer edge of the iron core, the mass of the silicon 
layer, the outer edge of the silicon layer, and the mass of the oxygen layer, respectively. The 
edge between each layer is defined as the radius where the most abundant element shifts to 
one another (see, Figured]). The mass of oxygen layer for the 15 M models of LC15 and 
WHW15 (Mo in the table) is larger than the other progenitors (i.e., WW15 and WHW11) 
and their oxygen layers (denoted by -Rsi/o) are positioned much closer to the center, so that 
they can touch to the supernova shock in a shorter timescale after bounce (before the neu- 
trino luminosity gets smaller with time). As one would anticipate, the impacts of nuclear 
burning are most remarkable for the LC15 progenitor as we will show in the later sections. 



Table 1: Summary of progenitor models and their composition features (see text for the 
definition of each quantity). 

Model M^ e R^> W Q 

(Me) (Me) (10 3 km) ( M Q ) (10 3 km) ( M ) 

LC15 13.4 1.44 1.31 0.221 2.22 0.814 

WHW15 12.6 1.55 1.96 0.124 2.97 0.943 

WW15 15.0 1.42 1.31 0.436 7.44 0.649 

WHW11 10.8 1.24 1.00 0.168 3.74 0.289 
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3. Results 




200 400 600 800 

Time [ms] 



Fig. 2. — Time evolution of the shock radii for model LC15 with different initial neutrino 
luminosities {L u q 52 in unit of 10 52 erg s _1 ) and the decay time (t^ in unit of s). The red 
and blue line corresponds to the results with and without the energy feedback from nuclear 
reactions, respectively 

In section 3.1, we start to investigate how the energy feedback from nuclear burning 
could affect the postbounce dynamics in ID simulations. Then we study how the nuclear- 
burning impacts are sensitive to the progenitor models, namely by the precollapse structures 
and their composition profiles (section 3.2). In section 3.3, we then move on to discuss how 
nuclear burning would affect the 2D dynamics. 

3.1. Impact of Nuclear Burning in ID simulations 

Relying on the light-bulb scheme in this study, the destiny of the stalling bounce shock 
(whether it will revive or not) depends simply on the two parameters; the input neutrino 
luminosity L u0 and the decay time td (see Equation ([6]))B Note in the following that we 
characterize models as (Lvo.52,^) = {x,y) for convenience, in which the luminosity and the 
decay time is x x 10 52 (erg/s) and y (s), respectively. Figure |2] shows comparisons of the 
postbounce shock evolution in ID LC15 models depending on the two parameters with (red 
lines) or without (blue lines) the a network calculation. 



4 Without seeing a shock revival in ~ 1 s postbounce, we call it as "non-exploding" in this study 
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For the chosen three cases in the figure; (L u0t52 , td) = (2.2, 5.0) (solid line), (2.0, 5.0) (dot- 
ted line), and (2.2, 2.0) (dashed line), all of the models with nuclear burning (red lines) show 
a shock expansion leading to explosions finally, while among the models without nuclear 
burning the solid blue line only exhibits a shock revival. This is the case of relatively higher 
luminosity with longer decay time ((-£^0,52, td) = (2.2, 5.0)). Note that in all the six cases (3 
different parameters x 2 (with/without burning)) in Figure 2, the bounce shock firstly stalls 
and then transits to a passive shock with no positive velocity (seen as a shrink of the shock 
radii at around ~ 50 ms postbounce in Figure 2). And only after that, the additional energy 
gain due to nuclear burning acts to bifurcate the path of the passive shock, namely whether 
the shock experiences recession afterward (for all the blue lines in Figure 2) or expansion 
(for red lines) with different revival timescales depending on the input neutrino parameters. 

As seen from Figure 2, larger input neutrino luminosity and shorter decay timescale 
(unsurprisingly) lead to more easier explosions. More importantly, by comparing dotted red 
with dotted blue line ((Lj,o,52, td) = (2.0, 5.0)), the shock is shown to shift from recession 
to expansion by the inclusion of nuclear burning. In the case of more luminous neutrino 
(-£-1/0,52 = 2.2), the trajectories of the shock are observed to be rather similar when the effect 
of nuclear burning is taken into account (compare solid red with dashed red line). 

In the following, we elaborate on how and why the shock expansion is affected by nuclear 
burning as observed in Figure 2. Figure [3] and H] show the mass-shell trajectory of models 
LC15 with the different parameter set of (L u05 2,td) = (2.2,2.0) and (L U Q 52 ,td) = (2.0,5.0), 
respectively. Note here that the former and latter case corresponds to the dashed and 
dotted line in Figure 2. Without nuclear burning (left panels in Figure [3] and H]) , the stalled 
shock, albeit showing several oscillations, does not turn into expansion (see also Figure [2]). 
With nuclear burning, by looking carefully at the right panels of Figures 3 and 4, the shock 
expansion can be seen to take place in the following two ways, whether the shock front passes 
through the Si-rich layer (see the behavior of the thick red line in the greenish region in the 
right panel of Figure [3]) or later it touches to the O-rich layer (e,g., the shock in the reddish 
region in the right panel of Figure H]). For the latter case, the bounce shock firstly stalls as 
in the non-burning model (compare the left with the right panel in Figure 4), but then the 
shock front says good-bye to the non-burning case when it encounters with the O-rich layer. 



200 400 600 800 
Time [ms] 



200 400 600 800 
Time [ms] 



Fig. 3. — Evolution of model LC15 with a parameter set of {L u0 ^2, td) — (2.2, 2.0) visualized 
by the mass-shell trajectories. The thick red line starting at t = denotes the position of 
the shock. Both cases either without (left) or with (right) the energy feedback from nuclear 
reactions are shown. The regions colored by gray, green, and red correspond to the iron, 
silicon, and oxygen layers, respectively. Thick gray lines correspond to the mass coordinates 
from 1.3 to 1.8 M & with every 0.1 M & (thin gray lines with every O.O2M ). 




Fig. 4. — Same as Figure 4 but for the parameter set of (L u0 ^ 2 ,td) = (2.0,5.0). 
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To look more in detail how the nuclear burning contributes to the shock acceleration, 
Figure 5 shows the radial velocity profiles and the composition distributions for model LC15 
with the parameter set of (£,,0,52,^) = (2.2,2.0) (the same parameter set as in Figure 3). 
At £ p b = 150 ms (top left panel), the shock front is in the silicon-rich layer (seen as a rapid 
decline of the green line in the mass fraction plot). Behind the shock front, heavier elements 
are synthesized as shown. The nuclear energy released by the silicon burning heats the 
material behind the shock, making the sign of the velocity slightly positive there (compare 
the velocity profiles with and without nuclear burning in the top left panel). The difference 
between the velocity profiles with versus without nuclear burning becomes outstanding when 
the oxygen-rich layer starts to touch the shock front (t p b > 250 ms). The left panel of Figure 
[6] shows the evolution of the explosion energy with (red line) /without (blue line) burning 
and also that of the net energy released via nuclear reactions (green line) . After the silicon 
burning starts to feed the energy behind the shock in addition to the neutrino heating (in 
the gain region, e.g., t p b = 150 ms, see the top left panel in Figure 5), the explosion energy 
deviates from the one without burning (compare red with blue line in the left panel in Figure 
6), which is also clearly visible in the shock evolution (compare red with blue line in the 
right panel in Figure 6). From the left panel of Figure 6, the total amount of 3.1 x 10 50 erg 
is shown to be released through nuclear burning in this case, lifting up the explosion energy 
to be 5.0 x 10 50 erg. 

As we already mentioned, the oxygen burning predominantly triggers the shock expan- 
sion for the parameter set taken for Figure 4. But also in this case, the silicon layer is 
shown to be burned as a heating source (top left panel of Figure 7), which is the reason 
that the shock position becomes larger compared to the non-burning model (see the right 
panel of Figure 8). When the shock front begins to swallow the oxygen layer at ~ 400 ms 
postbounce (the right panel of Figure 4), the fresh fuel supplies energy to assist the shock 
expansion (see, from top right, bottom left, to bottom right panels of Figure 7). If not for 
the energy gain, the stalled shock does not revive no earlier than t = 750 ms as seen from the 
right panel of Figure 8. Even with the aid of nuclear burning, the explosion for this model 
(left panel of Figure 8) is weaker (< 10 50 erg) comparing to the more luminous models (left 
panels of Figures 6 and 9). This suggests that nuclear burning has a secondary impact on 
the explosion mechanism in the sense that it can assist explosions only when the neutrino 
heating is working enough strong to push the weak shock to the fuel layers. 

In this section, we have elaborated on the features of the burning affects by taking the 
three different models as a reference. In the next section, we proceed to study a systematic 
trend by showing the parameter region spanned by L u0 and (as well as its progenitor 
dependence), where nuclear burning is able to affect the neutrino-driven explosions. 
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Fig. 5. — Snapshots of velocity profile (top) and composition distribution (bottom) for model 
LC15 with L u0 = 2.2 x 10 52 erg s _1 and t d = 2.0 s at selected postbounce epochs (£ p b = 
150, 250, 350, and 400 ms). Solid and dotted line in the top of each panel shows the velocity 
profile with or without nuclear burning, respectively. In the bottom part, distributions of 
representative elements are shown. Note that the abundances of neutron (n) and proton 
(p) in the NSE core (regions below the horizontal dashed line in the abundance plot) are 
estimated from Shen's EOS and the abundances of the other elements outside are calculated 
from the nuclear network calculation. 
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Fig. 6. — For the same model as Figure El the evolution of explosion energy (left) and the 
shock radius (right) are shown, respectively. In the left panel, the green line is the energy 
released by the nuclear burning and the explosion energy without nuclear burning (blue) is 
also shown. Here the explosion energy is defined as the sum of the kinetic, thermal, and 
gravitational energy of fluid elements with positive radial velocity. 
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Fig. 7. — Same as Figure 5 but for the parameter set of L u q = 2.0 x 10 52 erg/s and td = 5.0 
s. 
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Fig. 8. — Same as Figure 6 but for the parameter set of L„ 
s. 
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Fig. 9. — Same as Figure 6 but for the most energetic case, (L U Q 52 ,td) = (2.2,5.0), among 
the three examples shown in Figure 2. The explosion energy with nuclear burning is about 
8.0 x 10 50 ergs at t p b = 465 ms and still keeps rising. The net energy released via nuclear 
reactions at this time is ~ 3.0 x 10 50 erg, occupying a significant fraction (~ 40%) of the 
explosion energy. For the model without nuclear burning, the explosion energy is ~ 0.8 x 10 50 
erg at that time and closely saturates to be E cxp ~ 2.3 x 10 50 erg afterward. 
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3.2. Progenitor Dependence 
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Fig. 10. — Parameter maps of the initial neutrino luminosity (L u0 ) and its decay time (td) 
that separates the non-exploding regime (blueish region) from the exploding one (reddish 
regime) in ID simulations for the four different progenitors. A horizontal yellowish region 
in-between (clearly visible for the LC15 progenitor; top left) shows the parameter region in 
which ID explosions are obtained when the network calculation is performed. 

Figure [10] shows the parameter region on the (L^q, td) plane, in which ID explosions are 
obtained only if the network calculation is performed (the yellowish region sitting between the 
exploding and the non-exploding regime). The burning-mediated regime is clearly visible 
only for the LC15 progenitor. As already mentioned in section 2.2, this is because this 
model possesses a massive oxygen layer and the oxygen shell is positioned closest to the 
center among the progenitors taken in this study. 



The area of the yellowish region in Figure 10 is not so big even for the LC15 progenitor, 
which suggests again that nuclear burning has a secondary importance because it can work 
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Fig. 11. — Ratio of the explosion energy with nuclear burning to that without for model 
LC15 on the (L u0 ,td) plane. The dashed line is the boundary between the exploding and 
the "nuclear-burning- aided" regime in Figure [10J In this plot, the radius of each circle is 
taken to be proportional to the ratio of the explosion energy, where the ratio (the size of the 
circle) becomes smaller to be approaching to unity for models with higher input luminosity 
and with longer decay time (from lower left to upper right direction in this map). 



- 19 - 



only when the neutrino-driven explosion takes place very marginally. In the case of the 
marginal explosions (which may be often the case in recent first-principle CCSN simulations), 
however, it should be emphasized that the inclusion of nuclear burning does lead to the 
enhancement of the explosion energy. Figure 11 shows the ratio of the explosion energy 
with nuclear burning to that without for model LC15 on the (L u o,td) plane. In fact, the 
explosion energy is shown to be remarkably enhanced in the case of marginal explosions (i.e., 
the size of the circle becomes larger on the critical curves (dashed line)), and getting small 
for large L u and t<i, in which explosions are predominantly triggered by neutrino heating. 
As repeatedly mentioned so far, these features due to nuclear burning are only remarkable in 
the LC15 progenitor. Nevertheless, it is worth mentioning that the shock extent even for the 
WW95 progenitor (which is the most pessimistic case, see the bottom left panel of Figure 
10) becomes bigger for models with nuclear burning compared to those without (Figure [T2~j) . 

The critical luminosity for explosions can be read from the y-axis in Figure 10 in the limit 
of long td (namely, approaching to a constant neutrino luminosity), which corresponds to 2.7 
(WHW15), 2.0 (LC15), 1.9 (WW15), and 1.2 (WHW11) in unit of 10 52 erg/s, respectively. 
The critical luminosity becomes smallest for model WHW11 (mainly due to a compactness 
of the precollapse core and its tenuous envelope). When the input luminosity is taken below 
the critical curves, nuclear burning cannot alone drive explosions because the shock needs 
to be kicked out of the iron core firstly by neutrino heating (i.e., the shock revival due to 
neutrino heating is preconditioned to enjoy the assistance from nuclear burning). 
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Fig. 12. — Time evolution of the shock radii for several ID models employing the WW15 
progenitor. The shock moves farther out for models including nuclear burning (red lines) 
compared to those without (blue) 



-21 - 



3.3. 2D Results 

To see clearly the impacts of nuclear burning in our 2D simulations, we choose to employ 
the LC15 in the following. 

Figures [T3] and HH show entropy evolution (left- hand- side in the left panel) with the 
abundances of representative elements of oxygen, silicon, and nickel (from the right-hand- 
side in the left panel to the right-hand-side in the right panel) for two sets of neutrino 
parameters at selected postbounce epochs (£ p b = 100(top), 200(middle) and 300 (bottom) ms 
postbounce), respectively. Small- and large-scale inhomogeneities in the entropy plots come 
from neutrino-driven convection and the SA SI, both of which lead to more easier explosions 



in 2D than ID (e.g., Marek fc Jankal 120091 : lOhnishi et all 120061 : Murphy fe Burrows 



2008|) 



Reflecting the stochastic motions of the expanding shocks, the way how the (anisotropic) 
shock surfaces touch the nuclear fuel (in the shape of spherical shells) changes from models 
to models in 2D. In the case of Figure 13, the anisotropically expanding shock firstly reaches 
to the oxygen layer near in the vicinity of the north pole at t ~ 200 ms, which is seen as 
a penetration of the blueish region to the reddish region there (;oxygen layer (middle left)). 
Simultaneously, heavy elements like silicon and nickel are shown to be synthesized there 
(middle right), which helps to push the burning material preferentially along the direction 
for the moment (see the reddish region of the nickel abundance in the bottom right panel). 
In a more luminous case assuming larger luminosity and longer decay time (Figure 14), the 
SASI seemingly develops only weakly, barely deforming the the shock surface (essentially 
spherical), although neutrino-driven convection does make the postshock entropy distribution 
non-uniform. 

As shown from Figure [T5| nuclear burning does assist 2D explosions similar to ID, but 
the energy difference (here ~ 0.2 x 10 51 erg) is generally smaller in 2D than in ID (compare 
with Figures and [9]). The comparison of the energy gain due to nuclear burning between 
ID and 2D models is more clearly shown in Figure 16. In the figure, the difference of the 
explosion energy (AE ew ) with and without nuclear burning is shown as a function of the mass 
of the oxygen layer that is burned due to the shock passage (Mo,bm) both in ID (open circles) 
and 2D (closed circles) models for the LC15 progenitor. The slope of AE C ^ P as a function 
of M . b rri is shown to be milder in 2D than in ID. This may be because neutrino-driven 
convection and the SASI in 2D models (as indicated by entropy distributions in Figures 
[TBI and m]) could enhance the neutrino heating efficiency, which makes the contribution of 
nuclear burning relatively smaller compare to ID models. 

Finally Figure [17] shows a parameter map in 2D for the LC15 progenitor (compare 
Figure [TU] in ID). As expected, 2D hydrodynamics leads to more easier explosions compared 
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Fig. 13. — 2D distributions of entropy and representative elements (0, Si, and Ni) of model 
LC15 for (L u0 ,52,t d ) = (2.7,1.1) at t ph = 100 (top), 200 (middle), and 300 ms (bottom) 
postbounce, respectively. In this case, the anisotropically expanding shock firstly reaches 
to the oxygen layer near in the vicinity of the north pole at t ~ 200 ms, which is seen as 
an aspherical penetration of the blueish region to the oxygen layer (reddish region (left)). 
Looking at the corresponding region in the middle right panel, heavy elements like silicon 
and nickel are shown to be synthesized. 




Fig. 14. — Same as Figure [131 but of a 2D model with (L U Q 52 ,td) = (3.0, 1.1). In this case 
neutrino irradiation is very strong and its effective heating leads to nearly spherical shape. 
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Fig. 15. — Time evolution of explosion energy (red line) and net nuclear burning energy 
(green) for LC15 model with L v $ = 2.4 and t& = 3.0. Explosion energy of the case without 
nuclear burning (blue) is also shown. 
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Fig. 16. — Shown is the differences of explosion energy between models including and not 
including the effects of nuclear burning as functions of the consumed mass of oxygen. Re- 
sults from spherical (open circles) and 2-dimensional (filled) simulations for LC15 model are 
plotted. Lines are for linear fit of both data. 
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to ID (see the dotted lines which are the critical curves in ID). More importantly, the 
yellowish region does form in 2D models for the LC15 progenitor. It should be interesting 
to perform multi-D (radiation-hydro) simulations with nuclear network calculation for the 
virgin progenitor model. Given a number of uncertainties in the pre-supernova evolution 
calculations (such as in treatments of multi-D effects, mass- loss, and weak interactions), 
our results would indicate that the consideration of nuclear burning remains as one of the 
unignorable ingredients in the first-principle simulations as a potential remedy to cure these 
underpowered symptoms. 




Fig. 17. — Same as Figure [TU1 but for 2D simulations. The dashed lines represent the critical 
curves in ID (compare Figure fTTJT) . Note that we adopt coarse mesh points in the polar 
direction (32 uniform grids), so that we can perform 2D simulations for 174 models in total 
to make this parameter map. 



4. Conclusion 

We explored the potential impacts of nuclear burning on fostering the onset of the 
neutrino-driven explosions of core-collapse supernovae. By changing the neutrino luminosity 
and its decay time to obtain parametric explosions in ID and 2D models with or without a 
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13-isotope a network, we studied how the inclusion of nuclear burning co uld affect the post- 
bounc e dynamics for four progen itor models; three for 15.0 M stars of iLimongi fc Chieffi 
j2006h.lWooslev fc Weaver] Jl995h . and lWooslev et all ([SoJ, and one for an 11.2 M star of 
Woosley et al.l (120021 ) Our results showed that the energy gain due to nuclear burning of in- 
falling material behind the shock can energize the shock expansion especially for models that 
produce only marginal explosions in the absence of nuclear burning. These models enjoy the 
assistance from nuclear burning typically in the following two ways, whether the shock front 
passes through the silicon-rich layer, or later it touches to the oxygen-rich layer. Depending 
on the neutrino luminosity and its decay time, the explosion energy was found to increase 
up to a few times 10 50 erg for models with nuclear burning compared to the corresponding 
models without. The difference in the explosion energy becomes generally smaller in 2D than 
in ID, because neutrino-driven convection and the SASI in 2D models enhance the neutrino 
heating efficiency, making the contribution of nuclear burning relatively smaller compared 
to ID models. It was pointed out that these features are remarkable only for the LC15 
progenitor both in ID and 2D, which possesses a massive oxygen layer with its inner-edge 
radius being smallest among the employed progenitors, which makes the timescale shorter 
for the shock to encounter the rich fuel. Considering uncertainties in progenitor models, our 
results indicate that nuclear burning should still remain as one of the unignorable ingredients 
to foster the onset of neutrino-driven explosions. 

An important extension of this study would be going to 3D simulations with the same 
setting. Without the coordinate symmetry axis, the anisotropy of the expanding shock 
is expected to be much milder in 3D than in 2D. If this would be the case, the impacts 
of nuclear burning in 3D could be more important tha n in 2D, provided that the shock 
morphology in 3D would be more close to that in ID ( iHanke et al.ll201ll ). It should be 
interesting to employ the LC15 model, which has been unexplored by 2D and 3D radiation- 
hydro simulations including a network calculation. Suppose that nuclear burning impacts 
the explosion dy namics and its blast morphology, their consequences on the multi- messenger 
observables (see, Kotake et al.l I2012L for a recent review) su ch as kick/spin of pulsars (e.g ., 



Wongwathanarat et al.ll2010h and nucleosynthetic yields (e.g.. lArcones fc Thielemann! (120121); 



2009 



Fuiimoto et al.l (120111) ) as well as neutrino and gravitational-wave signals (e.g 



Ott et al. 20121: Lund et al. 2010: Muller et al. 2012: Kotake et al. 2009a.b 



Marek et al. 



2011 



2012 



might deserve further investigation. We hope that our results would give a chance/momentum 
for CCSN researchers to reexamine the nuclear-burning effects on their own results. 
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